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Abstract 


In this study, we apply a two-dimensional variational analysis method (2d-VAR) to select 
a wind solution from NASA Scatterometer (NSCAT) ambiguous winds. 2d-VAR determines 
a “best” gridded surface wind analysis by minimizing a cost function. The cost function 
measures the misfit to the observations, the background, and the filtering and dynamical 
constraints. The ambiguity closest in direction to the minimizing analysis is selected. 

2d-VAR method, sensitivity and numerical behavior are described. 2d-VAR is compared 
to statistical interpolation (01) by examining the response of both systems to a single ship 
observation and to a swath of unique scatterometer winds. 2d-VAR is used with both NSCAT 
ambiguities and NSCAT backscatter values. Results are roughly comparable. When the 
background field is poor, 2d-VAR ambiguity removal often selects low probability ambiguities. 
To avoid this behavior, an initial 2d-VAR analysis, using only the two most likely ambiguities, 
provides the first guess for an analysis using all the ambiguities or the backscatter data. 

2d-VAR and median filter selected ambiguities usually agree. Both methods require 
horizontal consistency, so disagreements occur in clumps, or as linear features. In these 
cases, 2d-VAR ambiguities are often more meteorologically reasonable and more consistent 
with satellite imagery. 


1 Introduction 


Scatterometers are active radars designed to measure the backscatter or normalized radar 
cross section (NRCS) from the earth’s surface at moderate incidence angles. The advantage 
of using moderate incidence angles (20 — 70°) is that the major mechanism for scattering 
is then Bragg scattering from centimeter-scale waves, which are, in most conditions, in 
equilibrium with the local wind. Although the scatterometer winds are usually provided as 
neutral winds at some reference height, the measurement is physically most closely connected 
with surface stress (Brown 1986 [7]). 

Scatterometers have been mounted on stationary platforms, aircraft and satellites. Since 
1991, satellite-borne scatterometers have provided a wealth of wind data over the world 
ocean. Previous satellite scatterometers include Seasat-A Satellite Scatterometer (SASS) in 
1979; the active microwave instruments (AMI) on the first and second European Remote 
Sensing satellites (ERS-1 and ERS-2), during 1991-2000; and the NASA Scatterometer 
(NSC AT) on the first Advanced Earth Observing Satellite (adeos-1) during 1996-1997. 
Currently and since 1999, SeaWinds on QuikSCAT has been producing useful data. 

Scatterometers measure wind speed very accurately; Freilich and Dunbar (1999 [15]) 
report that NSCAT measured wind speed to within 1.3 ms -1 in an rms sense, while Stoffelen 
and Anderson (1997 [55]) suggest that the ERS-1 measurement error standard deviation is 
actually only 0.5 ms -1 . However, the wind direction from scatterometers is not uniquely 
determined. Wind speed and direction are inferred from a number of closely colocated 
(both temporally and spatially) radar measurements from a number of different azimuth 
angles. The measured NRCS, denoted a 0 , varies with the relative azimuth angle between 
the antenna and wind direction. An “upwind-crosswind” variation exists because small, 
wind-generated gravity waves on the ocean surface reflect more of the radar signal when 
wave crests are perpendicular (i.e., “upwind”) to the radar antenna than when wave crests 
are parallel (i.e., “crosswind”). The “upwind-downwind” signal is very small, because when 
fetch is not limited, the orientation and amplitude of these waves is nearly the same if the 
wind direction is shifted by 180°. 


1 



Typically several a 0 measurements are grouped into 25 or 50 km wind vector cells 
(WVCs) for wind retrieval. Modeled a 0 values are computed from a geophysical model 
function which requires spacecraft position and orientation (known) and wind speed and 
direction (unknowns) as inputs. The normalized squared differences between measured and 
modeled <r° values are summed as part of the objective function of a maximum likelihood 
estimator (MLE). Then, wind speed and direction are varied to minimize the difference be- 
tween measured and modeled <r° values. However, over the range of possible wind speeds and 
directions (i.e., 0-50 ms’ 1 and 0 - 360°), the MLE produces from two to six local maxima, 
and it is not possible to determine the wind uniquely from the MLE alone. The wind vectors 
corresponding to these maxima are called “ambiguities”. For operational numerical weather 
prediction (NWP), and for other purposes it is necessary to choose one ambiguity at each 
location in a meteorologically consistent manner. Once we have chosen a single ambiguity, 
we call the result a “unique” wind. 

Median filter techniques have been used for ambiguity removal, but may fail at swath 
edges and in regions of strongly curved or sheared flow. For NSC AT, the Jet Propulsion Lab- 
oratory (JPL) used a median filter initialized by the ambiguity with the highest MLE. This 
method was completely objective and autonomous. Later, better performance was obtained 
by initializing the scatterometer wind field with the ambiguity closest to the operational 
gridded surface wind field analysis from the National Centers for Environmental Prediction 
(NCEP, formerly the National Meteorological Center or NMC). This product is referred to 
as the nudged wind product since an outside source is used to initialize (or “nudge ) the me- 
dian filter. This method is also completely objective but depends on information in addition 
to the scatterometer measurements. The performance of the median filter method has been 
quite successful (Gonzales and Long 1999 [18]) but is not based on any physical model of 
wind flow. Consequently, the median filter occasionally selects physically unrealistic winds. 
Patches of unrealistic winds in the scatterometer data might have a disastrous impact on the 
quality of analyses and NWP forecasts. Careful quality control (QC) is required, of course, 
but is problematic if the model winds were used in the ambiguity removal. 
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This paper reports on the use of a two-dimensional variational analysis method (2d-VAR) 
for ambiguity removal. 2d-VAR generates a gridded surface wind analysis by minimizing an 
objective function, which is a weighted sum of several constraints on the difference between 
the analysis and the background, and between the analysis and the data. While the “nudged’ 
median filter uses a priori information for initialization, 2d-VAR uses such information both 
for initialization and as part of the background constraints. 2d-VAR was originally described 
by Hoffman (1982 [22]) to analyze Seasat scatterometer (SASS) data. In this first formu- 
lation, only observations from conventional platforms (ships, buoys, radiosondes) and from 
SASS were combined with an a priori or background wind field. 2d-VAR was extended to 
include smoothness (i.e. filtering) constraints and a dynamical constraint by Hoffman (1984, 
hereafter H84 [23]). 

Since H84 2d-VAR has undergone several developments. These include: 

• Optional use of 12 point bicubic interpolation to the “obs” locations (§B). 

• Use of <7° values (§6.a.l). 

• Reformulation of the constraint terms (§3.b). 

2d-VAR may combine both wind retrieval and ambiguity removal by analyzing the <7° data 
directly. For both ERS-1,2 and NSC AT, in addition to the possibility of analyzing a 0 data, 
the ambiguities may be treated in the same manner as SASS data were treated by H84, 
or unique winds may be treated as buoy “obs” . 2d-VAR has been used to produce several 
long running data sets for oceanographic research. First, 2d-VAR was adapted for use with 
time binned wind observations by FSU researchers including Legler et al. (1989 [33]), Meyers 
et al. (1994 [38]), Jones et al. (1995 [28]), and Pegion et al. (2000 [45]). Later, 2d-VAR was 
extended to SSM/I wind magnitude data by Atlas et al. (1991, 1996 [2, 3]). We will review 
the 2d-VAR formulation in section § 3, and extensions for scatterometer data in § 6. 

2d-VAR, in common with other 3d-VAR and 4d-VAR systems, is based on a Bayesian 
viewpoint. In fact the 2d-VAR cost function may be written as J = J b + J a , where 

• J is the total cost function. 
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J b is the background cost function. 


• J 0 is the “obs” cost function. 

Since J is a function of the surface wind field, minimizing J produces a gridded surface wind 
analysis that simultaneously closely fits the data and background (H84). As in any analysis 
method, specification of the errors of the data and the background are vitally important. 
Our implementation of 2d-VAR employs a heuristic model of forecast error statistics. J b is 
defined below in (3) as the weighted sum of several constraints. Each constraint function, 
J m , is a scalar which measures the the difference between the analysis and background 
using the mth constraint. The constraints include filtering and dynamical consistency. The 
combination of these constraints mimicks the usual J b term in 3d-VAR or 4d-VAR (e.g. 
Thepaut et al. 1993 [57]). Therefore 2d-VAR as described here may be considered to be 
a specialization of 3d-VAR to ocean surface wind. Wahba and Wendelberger (1980 [58]) 
showed that under certain conditions 3d-VAR and 01 are the same. Therefore 2d-VAR is 
closely related to 01 as well. Although the formulation of 2d-VAR described here is very 
different from standard 01 or 3d-VAR, analysis results can be made similar by appropriate 
choice of lambda weights (A,, multiplying factors which control how much each function 
contributes to the total objective function). (See §5.) Thus 2d-VAR is a form of 01 in 
which the correlations are parameterized with a small number of degrees of freedom. After 
subjective tuning, the 2d-VAR heuristic model of forecast error statistics produces results 
similar to conventional error models, but has the potential advantage of being specified by 
a small number of parameters. In principle, these parameters might vary with geographic 
location, with synoptic situation and might even be tuned online, using current or recent 
data. As a result, 2d-VAR is an ideal candidate for on-line estimation of parameters as Dee 
et al. (1998 [11, 12]) have proposed. This feature of 2d-VAR has been exploited by Pegion 
et al. (2000 [45]) to objectively tune these parameters. In the work described here, these 
parameters are chosen subjectively based on experimentation. However some of our results 
address how r to tune these parameters to get similar results when changing the analysis 
resolution, data density, or data representation (see § 7.) 
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The outline of the rest of this paper is the following. First we define the NSC AT data and 
the nominal processing of the NSCAT data (§2). We then describe 2d-VAR for conventional 
data (§3). Experiments with a single ship observation then illustrate the numerical behavior 
and the sensitivity of 2d-VAR (§4). These experiments demonstrate that the subjectively 
tuned 2d-VAR approximates 01. Idealized 01 experiments with a swath of unique NSCAT 
winds are presented to demonstrate the close relationship between 2d-VAR and 01 (§ 5). The 
use of 2d-VAR for ambiguity removal for NSCAT is then described (§ 6). Experiments using 
ambiguities and cr° values are compared. With proper tuning the results are similar (§7). 
We have improved the quality of our ambiguity removal by using a first stage in which only 
“dual-ambiguities” are processed (§8). Examples of our dual-QC method in 2d-VAR are 
shown in §8. a. Then §9 is a conclusion. In a companion article (Henderson et a/.2001 [21]) 
the dual-QC 2d-VAR. is applied to the entire NSCAT mission and is tested as an alternative 
to the standard ambiguity removal schemes in impact experiments using the GEOS-2 data 
assimilation system. 

2 NSCAT data 

The focus of this study are winds retrieved from NSCAT. NSCAT was launched aboard the 
ADEOS-1 spacecraft on 16 August 1996 from Tanegashima Space Center in Japan. Unfortu- 
nately, after only nine months of operation, the ADEOS-1 solar power array failed. ADEOS-1 
had a nearly polar Sun-synchronous orbit at a height of ~ 800 km, with a period of ~ 100 
minutes. 

NSCAT made observations covering two 600 km swaths, one on either side of the space- 
craft, separated by a gap of ~ 330 km corresponding to low incidence angles. On each side 
of the spacecraft, NSCAT operated three 3 m long, stick-like antennas at 13.995 GHz. The 
foreward and aft antennas were vertically polarized, while the mid antenna was both ver- 
tically and horizontally polarized. Thus NSCAT made observations in sequence from eight 
beams, i.e. from eight combinations of antennas and polarizations. While the foreward 
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and aft antennas are separated by 90° as in the ERS 1 design, the angles between the mid 
antenna and the others are 20 and 70°, whereas in the ERS— 1 design these angles are 45 . 
Complete descriptions of the NSCAT antenna subsystems, resolution, registration of ob- 
servations, measurement sequencing, hardware, and ground system are provided by Naderi 
et al. (1991 [39]). 

In the ground processing, the individual backscatter observations associated with the 
different beams are organized by a grid oriented along the satellite track. Nominally there 
are 4 or 16 backscatter observations in a 25 or 50 km WVC. The WVC coordinates along, 
and perpendicular to, the satellite track are called row and cell, respectively. Each cell is 
associated with a nearly constant incidence angle for each antenna, with the values being 
equal for the fore and aft antennas. Wind ambiguities are retrieved for each WVC contain- 
ing a sufficient number of quality controlled backscatter observations. The retrieved winds 
maximize a likelihood function locally. The likelihood function in turn depends on the geo- 
physical model function to compare a candidate wind to the backscatter observations (Chi 
and Li 1988, Offiler 1994 [8, 44]). Examples of the likelihood function are described in § 6.a.l. 

Geophysical model functions are empirical relationships used to relate the backscatter to 
the geophysical parameters, and are derived from colocated observations (Jones et al. 1977 
[30]). In current model functions, the backscatter depends nonlinearly on wind speed and 
direction. Models developed from aircraft flights have been used as baseline model functions 
for satellite instruments. However substantial post-launch modifications and refinements 
were needed for SASS (Boggs 1981 [5]) and for ERS-1 (Offiler 1994, Stoffelen and Anderson 
1997 [44, 55]). The prelaunch NSCAT model function, denoted NSCAT-0 or SASS2, w^as 
developed by Wentz et al. (1984 [60]) from the Seasat data set by matching the statistics 
of the observed a 0 and those simulated assuming that the ocean surface winds follow a one 
parameter Weibull distribution (Conradsen et al. 1984 [9]). The fitting procedure used six 
parameters at each polarization and incidence angle. The first postlaunch NSCAT model 
function, denoted NSCAT— 1, was developed using a combination of SSM/I, NCEP and 
European Centre for Medium- Range Weather Forecasts (ECMWF) winds colocated with 
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NSCAT cr° observations. A further refinement, denoted NSCAT-2, has been developed by 
Wentz and Smith (1999 [61]), based on colocations with SSM/I and ECMWF winds. 

There have been many theoretical attempts to derive model functions for scatterometers 
but they have not been as successful as the purely empirical model functions. A moments 
consideration of the complexity of the ocean surface under moderate to high wind conditions 
suggests that a purely theoretical approach is a difficult and ambitious undertaking. Recent 
examples include Janssen (1994 [27]) and Romeiser et al. (1997 [49, 50]). 

Note that the model function relates wind speed and direction to the backscatter. In this 
calculation the geometry of the observation is taken as a given. To retrieve scatterometer 
wind information, several backscatter measurements at the same earth location, but with 
different viewing geometries, must be used. The winds are obtained by optimally fitting 
these data. 

Because of the nonlinearity of the model function and the various sources of error, several 
wind vectors consistent with the backscatter observations are usually found (Price, 1976 
[47]). (See § 6. a.) These multiple wind vectors are called aliases in the early literature and 
are now generally referred to as ambiguities. With two beams, SASS typically retrieved four 
ambiguities. With additional beams, NSCAT and ERS often retrieve just two ambiguities 
with directions approximately 180° out of phase. The ambiguities are ordered by likelihood. 
The first ambiguity is the ambiguity with the highest likelihood; it is most consistent with 
the backscatter data. The second ambiguity is usually of nearly equal speed and consistency, 
but of nearly opposite direction. 

Since the wind retrieval is performed for each WVC individually and since there is little 
to distinguish the first two ambiguities, a wind field of first ambiguities will generally be very 
irregular with many wind reversals. Spatial filtering can be used to produce a horizontally 
consistent wind field and resolve the directional ambiguity. A variety of approaches have 
been used. These include the work of Hoffman (1984 [23]), Schroeder et al. (1985 [51]), 
Schultz (1990 [52]), Shaffer et al. (1991 [53]), Wentz (1991 [59]), Badran et al. (1991 [4]), 
Offiler (1992 [43]), Thepaut et al. (1993 [56]), Long (1993 [35]), Stoffelen and Anderson 
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(1997 [54]), Jones et al. (1999 [29]), Figa and Stoffelen (2000 [13]), and de Vries and Stof- 
felen (2000 [10]). Simulation studies for NSCAT established that a median filter operating 
autonomously, i.e. with no a priori information, would provide excellent ambiguity selection 
(Schultz 1991, Shaffer et al. 1991 [52, 53]). The median filter must be initialized or seeded 
with some initial choice of ambiguities. In autonomous mode the seed is the first ambiguity. 
In practice this choice did not perform as well as expected. For the nudged median filter, the 
seed is the ambiguity closest to a background obtained from an operational analysis or short 
term forecast. No further reference to the background is made. At the end of the process 
the chosen winds are consistent with backscatter observations in the WV C and horizontally 
consistent with neighboring WVCs, but are not necessarily close to the background. Dif- 
ferences between autonomous and nudged results show that the median filter has multiple 
solutions. Thus spatial filtering does not resolve all ambiguity in the sense that multiple 

horizontally consistent wind fields are plausible. 

It should be noted that there are several NSCAT data products available. These differ in 
resolution (25 or 50 km), geophysical model function (NSCAT-0, NSCAT-1, or NSCAT-2), 
and ambiguity selection procedure (autonomous or nudged median filter). All JPL NSCAT 
data sets use the median filter described by Shaffer et al. (1991 [53]). The 25 km resolution 
data sets based on the NSCAT-2 model function and the nudged median filter are the most 
current. Some of the preliminary results presented here however use other versions of the 

NSCAT data. 

3 2d-VAR method 

2d-VAR finds the minimum of an objective function. The objective function depends on 
the surface wind field. Given any wind field the objective function measures the difference 
between the given field and the observations, and between the given field and the a priori or 
background wind field. The minimizer of the objective function is therefore the wind field 
which optimally fits the observed and a priori data. 
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The order of the steps taken in 2d-\'AR is flexible and the subject of ongoing tuning. The 
nominal and minimal setup is to: 

1. Define the geometry of the analysis grid, and on this grid define the background wind 
field. 

2. Set the inital estimate for the analysis equal to the background. 

3. Read the “obs” data sets, making necessary conversions and performing gross quality 
control (QC). 

4. QC the “obs” based on simulated values from the background. 

5. Minimize the objective function J . 

6. Write the analysis wind field and simulated “obs” values from the analysis. 

This basic operation may be modified in several ways. 

• The initial guess may be an arbitrary wind field. 

• After a preliminary analysis, the analysis and background can be interpolated to a 
finer grid. The minimization then restarts, using the newly interpolated analysis as 
the initial estimate of the fine scale analysis. 

• After a preliminary analysis, all data may be subjected to a second QC. This accepts 
some of the data which were initially rejected but are corroborated by near-by data. 
The minimization might then restart. 

In §8 a small modification of this basic setup provides the dual QC approach. 

In order to minimize the objective function J we must be able to calculate it for any 
physically reasonable value of the control variable vector X , which here represents the sur- 
face wind field. In the current formulation X contains those values of the gridded wind 
components (u, v) which are independent and wdiich are free to vary. Some points in the 
gridded wind field are not independent. In a global field points at longitudes 0 and 360 
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are the same. Not all points in the gridded wind fields are allowed to vary. If there are fixed 
points, they are not part of the control vector. The fixed points are nominally set equal to 
the background, but in general they are arbitrary. In particular for multiprocessing using 
domain decomposition the fixed points may be obtained from another processor’s previously 
computed solution. A detailed description of the gridded wind field representation is given 
in § A. 

Efficient minimization also requires the calculation of the gradient of J, i.e. dJjdX. 
This is accomplished by the adjoint of the routines which calculate J (Hoffman et al. 1992 
[26]). The minimization procedure we use is a version of the limited-memory quasi-Newton 
algorithm (Liu and Nocedal 1989 [34]), which was described and evaluated by Navon and 
Legler (1987 [40]). The stopping criteria is either 

|G| < emax(l, |X|), (1) 

or, until Nj function evaluations, whichever occurs first. Here 

• G is the gradient of J with respect to X. 

• e is the gradient test stopping parameter. 

• Nj is the function call stopping parameter. 

We have been unsuccessful at preconditioning the minimization. Typically several hun- 
dred iterations are required to satisfy our very strict convergence criterion. In a typical 
experiment the rms difference in the u or v wind components between the solution after 250 
and after 225 function calls is ~ 0.02 ms -1 , with maximum differences of ~ 0.5 ms -1 . The 
reason for the slow convergence is due to the step by step expansion of the region of influence 
of the data (see Fig. 5). A regridding process speeds convergence (§4.b.2). 

The J used in this study may be written in the form, 

J — Jb + Jo: (^) 

where 

Jfj AvwM J VWM 4” A l ap J LAP 4~ A j-, j y J D I V 4~ AyQpt/yOR 4” Adyjsj *7dYN T (3) 
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^conWconv 


+ ^AMB^AMB + A SPD J S PD + A 


Here, the weights, denoted A ,, control the amount of influence each constraint has on the 
final analysis. As mentioned in the Introduction, the A* are chosen subjectively in this study, 
but techniques to tune them objectively have been developed. Values of the A* which were 
used for various experiments are reported below. It should be noted that these are weak 
constraints and therefore the Aj are not Lagrangian multipliers. Table 1 contains a summary 
of the constraints used. Details of the formulation follow, or are given in the Appendices. 


3. a “Obs” functions 

In 2d-VAR, the misfit between data and analysis is measured by a objective function de- 
fined as a sum of squared errors over all data locations. Multiplying factors, called lambda 
weights and denoted by Aj, control how much each function contributes to the total objective 
function. For example, we might divide ships into special research ships and all others and 
buoys into Tropical Atmosphere-Ocean (TAO) buoys and all others, or we might lump all 
these data into the conventional platform class. Similarly we might divide the SSM/I data 
depending on the DMSP spacecraft designation (e.g., F8), or keep them all together. We 
expect the minimizing analysis to fit the data more closely for a given data type as the A 
weight for that data type is increased. 

In general, the total “obs” function, J 0 is the weighted sum of individual “obs” operators 
for each data class, but here we report experiments with only one type of data and some of 
the 

A 0 = (A C onvj A amb , A spd , A nrcs ) (5) 


are set to zero. Usually the “obs” function is given by the sum of squared differences between 
simulated (i.e., calculated from the analysis) and data values normalized by the estimated 
“obs” (and representativeness) error(s). For example, 

( u a ~~ Wo ) 2 . (w a ~ W 0 ) 2 

-'conv — „ 2 + c2 ' V J 
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Here 


• (ua,^) is the analyzed wind interpolated to the “obs” locations as described in §B. 

• ( u o ,v 0 ) is the observed wind. 

• s u is the wdnd component error standard deviation. 

If “obs” errors within the data class are correlated in a known way , these correlations can 
be accounted for as described by Thepaut et al. (1993 [56]), by setting J c onv — Z O Z. 
Here 

• Z is the vector of normalized departures. Z would contain an entry for each wind 
component of each “obs”. 

• O is the matrix of observation error correlations. 

In the 2d-VAR experiments reported here, “obs” error correlations are ignored. Currently 
s u is taken to be a constant 1 ms -1 and knowledge of the “obs” error standard deviation is 
absorbed into A C qnv- The special “obs” functions for scatterometer data are described later 
in § 6. a. 

3.b Background constraints 

As given in (3), the background function, Jb is a weighted sum of several constraints, each 
providing a measure of the difference between the analysis and background wind field. The 
simplest constraint measures the square of the vector wind magnitude: 

X»m = ~ ( ( V, - V t ) 2 dA. (7) 

d\'WM J A 

The effect of J vw m is to constrain the size of the analysis increments. 

The scaling parameter S V wm makes J v WM nondimensional and of convenient magnitude. 
[If these were not concerns 5 VW m might be absorbed into A VW m since the adjustable parameter 
A vwm multiplies J vw m in (3).] We define 5 V WM in terms of typical velocity and length scales, 

S VWM = V 2 L 2 = L 4 /T 2 . (8) 
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Here 


• L is the length scale [10 6 m — 10 3 km]. 

• V is the velocity scale [10 ms -1 ]. 

And from these we define a time scale for later use: 

• T is the time scale [L/V = 10 5 s]. 

The other background constraints which are summarized in Table 1 are given by 


Jlap = T 2 J A \V 2 {u a - u b )] 2 + [V 2 (u a - v b )] 2 dA, 

(9) 

Jmv = S / 4 [v • (k - 

(10) 

Jv OR = ~ /jv • k X (V a - V b )] 2 dA, and 

(11) 


(12) 


Here 

• (u b ,v b ) is the background wind. 

• £ is the relative vorticity. 

• k is the unit vertical vector. 

Details of the calculation of J b are given in the appendix (§ C). 

4 Response to a single ship observation 

Experiments with a single conventional wind observation are presented in this section. These 
experiments are similar to those in H84. However results of the current experiments are 
somewhat at odds with the earlier conclusions because H84 was overly optimistic about the 
convergence of the 2d-VAR minimization and did not always use sufficient iterations in high 
resolution experiments. Here we show that 
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• The 2d-VAR solution (after sufficient iterations) is only weakly dependent on grid 
resolution. (Fig. 3.) 

• The 2d-VAR solutions (each of which is the end of a convergent iteration) converge as 
resolution increases. (Fig. 4.) 

• The convergence rate of the 2d-VAR minimization as presently implemented is slow, 
especially as resolution increases. (Fig. 5.) 

• For practical purposes, a high resolution 2d-VAR analyses should be obtained as the 
last in a series of 2d-VAR analyses at increasing resolution, each initialized by the 
previous one. 

• At high resolution, the effect of interpolation errors becomes noticeable and bicubic 
interpolation should be used. (Fig. 9.) 

• The sum of two single “obs” solutions is qualitatively the same as the solution for the 
two “obs” together. (Fig. 10.) 

4.a Methodology 

Most of the results presented in this section are based on 2d-VAR analyses for a single ship 
at 42N, 50 W, observing 30 ms~ l , from 210°. For these experiments the background field is 
taken to be totally calm. The resulting analyses are essentially Green s functions, i.e., the 
response of the analysis to a single impulse. 

All analyses are done on a grid large enough that boundary effects should be small. The 
grid runs from 290E to 330E, and from 24N to 60N. For graphical output we thinned (or 
interpolated) the analysis grid to one degree resolution in the window 302E to 322E and 35N 
to 50N. 

Most results are for 

A = (AcoNVj AywM ; A LA p, A di v, AyOR, A DyN ) (80, 1, 4, 16, 4, 64) , 
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and the single ship “obs”. These A values are equivalent to the nominal case we have used 
before except that J v wm is given 1/4 of the weight used by H84. 

4.b Results 

Results for the current nominal case are given in Fig. 1. The response is nearly 20 ms at the 
observation location. The response is elongated in the direction of the wind observation, and 
at sufficient distance in the perpendicular direction, there is a return flow, most significantly 
to the northwest. 

The qualitative response for the nominal case is similar the response of ordinary statistical 
(or optimal) interpolation to a single wind observation. Fig. 2 show's the optimal interpolation 
result for the same situation, assuming the w'ind component “obs” standard deviation for 
ships is 2.2 ms~ l and using NOGAPS forecast error statistics (Goerss and Phoebus 1993 
[17]). Compared to the nominal case, these results are more symmetric, with the influence of 
the data extended further in the transverse direction and less in the longitudinal direction. 

4.b.l Effect of grid resolution 

Runs were made w'ith grid increments of 2, 1, 1/2, 1/4 and 1/8 using various methods of 
initialization and stopping criteria. 

For resolutions from 1 to 1/8° (Fig. 3 show's the results for 2 to 1/4 ) the solutions are 
nearly the same. That is, if everything else is fixed (weights, data, domain, ...) 2d-VAR 
produces the same answer independent of resolution, if it is allow'ed to iterate sufficiently. 
The 2° solution (Fig. 3) has somew'hat jagged contours and exhibits less return flow', but is 
similar to the higher resolution solutions. 

The implication of these findings is that for constraints like those we have been using, 
anything finer than a 1° resolution is unnecessary and only relatively small errors are in- 
troduced by using 2° resolution. Using a higher resolution 2d-VAR grid would not produce 
a higher resolution analysis, only more grid points. Weaker background constraints should 
produce an analysis more closely fitting high resolution data. 
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In Fig. 4 differences with respect to the 1/8° solution are displayed. Clearly, the solutions 
obtained are approaching a limit as the resolution is increased. Differences are concentrated 
at the “obs” location. These differences are small and in obtaining these results it is criti- 
cal that the 2d-VAR minimization has converged. Doubling the resolution, quadruples the 
degrees of freedom, quadruples the computation effort per iteration and quadruples the num- 
ber of iterations required to achieve comparable accuracy, thereby increasing the memory 
requirement by a factor of 4 and the timing requirement by a factor of 16. At high resolution, 
convergence to the true solution is quite slow, and for the results of Fig. 3 we used many iter- 
ations and computed the solutions in sequence of increasing resolution, using the previously 
computed half resolution solution as the initial estimate in each case. The iteration stopped 
according to (1) with c = 0.0001 and iVy = 1000. This occurred for resolutions of 2, 1, and 
1/2° after 73, 210, and 846 iterations, respectively. Higher resolution solutions terminated 
after 1000 iterations. 

4.b.2 Convergence of the minimization 

The reason for the slow convergence of the 2d-VAR minimization at high resolution is illus- 
trated by Fig. 5. This shows the results at 1/2° resolution, beginning with a zero initial 
estimate, after 25, 50, 100, and 200 iterations. 2d-VAR initially attempts to fit the ship 
“obs” locally. At any iteration, if a grid point currently has a zero analysis wind increment 
(relative to the background) and if all its neighbors do also, then all the constraints are 
perfectly satisfied and there will be no change during the next iteration. As a result the 
pattern of the solution slowly expands and grows during the iteration. 

Many iterations were needed to show r that the numerical solution is approaching a limiting 
(presumably true) solution. However for practical purposes it is sufficient to iterate 25 times 
at 2°, interpolate to 1°, iterate 25 more times, and continuing this process to the desired 
resolution. We call this process regridding. Regridding with only 25 iterates per resolution 
step provides agreement to the most accurate solution to within 1 ms 1 at resolutions of 
1/2° and better, and agreement to the converged solution for the same grid increment to 
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1 ms -1 or better at all resolutions. The latter result is illustrated by Fig. 6, and the former 
may be inferred from Fig. 6 and Fig. 4. 

4.b.3 Sensitivity to the A weights 

Analyses with varying A weights were made using 1/2° resolution. In these analyses (and 
in analyses presented in following sections), up to 1000 iterations were allowed after initial 
estimates were developed using the regridding process described in the previous section. 
In Table 2 the total number of function evaluations used is listed in the column headed 
Nf. This number is approximately 60 more than the number of iterations used in the final 
minimization. 

Fig. 7 shows selected results obtained by varying some of the A weights by a factor of 
4. The actual A weights used are listed in Table 2, except that in all cases A VW m = 1? arK l 
the background weight is effectively altered by changing all other weights in the opposite 
sense. Changing the “obs” weight (panel ‘ 1 obs ’ ’ +) changes the amplitude of the response, 
but does not change the shape or scale of the response. In fact the solutions for different 
A C onv, with everything else held constant, are remarkably well correlated in the sense that 
one solution is nearly equal to another solution times a constant. Increasing the filtering 
and dynamic constraints (panel Constraints +) gives a solution of larger spatial scale, with 
smaller amplitude and less return flow to the right and left of the wind “obs” . Decreasing the 
weight given the fit to the background (panel Background -) provides a generally bigger 
response of the same shape, but larger spatial scale, relative to the nominal case. 

The half-width at half-height of the 2d-VAR response to a single “obs” might be used to 
define the effective resolution of the analysis. Fig. 7 shows that this resolution is controlled by 
the weight given the filtering and dynamic constraints relative to the background constraint. 
On the other hand the amplitude of the 2d-VAR response can be controlled by varying A CONV - 
For the nominal case changing some of the A m by a factor of 4 results in a 25—30% change 
in amplitude or scale of the response (Table 2 and Fig. 7). 

Fig. 8 shows selected results obtained by setting some of the weights to zero. In general, 
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eliminating a constraint increases the amplitude of the analysis. Eliminating the divergence 
and vorticity penalty functions (panel w/o DIV, VOR) gives a solution of similar spatial 
scale, with larger amplitude and greater return flow to the right and left of the wind "obs” . 
Eliminating the dynamic constraint (panel w/ o DYN) results in a less elongated solution, less 
return flow and larger amplitude at the “obs” location. Eliminating both the divergence 
and vorticity penalty functions and the dynamic constraint (panel LAP only) gives a very 
symmetric pattern with no variation in wind direction. In this case the problem for each 
velocity component is effectively decoupled. In the cases without the dynamic constraint 
(the right column of Fig. 8), there is no dependence on the Coriolis parameter or indeed 
on location, and the solutions are symmetric about the line defined by the observed wind 
vector. 

4.b.4 Effect of interpolation 

In general gridded values of wind components are interpolated to the “obs' 1 locations. In 
the experiments described before, the “obs” location is coincident with a grid point. In such 
cases the interpolation weight is one for the coincident grid point, and zero for all other grid 
points. Effectively, no interpolation takes place in such cases. 

To examine the effect of interpolation we shifted the “obs” location to 42.25N, 49.75W. 
Analyses were created for varying resolutions and for bilinear and bicubic interpolation. 
Fig. 9 shows differences with respect to the solution for 1/4° resolution. Such differences 
show the combined effects of interpolation and resolution, and should be compared with 
Fig. 4. Note that for resolutions of 1/4° and finer, the “obs” location is coincident with a 
grid point, and the solution does not depend on the interpolation procedure used. The 1/4° 
solution is identical, but shifted by 1/4° north and east with respect to the nominal case. 

The linear interpolation produces a solution with larger analyzed values. This is antici- 
pated since bilinear interpolation on a grid cell, unlike higher order interpolations, can not 
produce any values outside the range of the gridded values surrounding the grid cell. 

In general the interpolation errors are the same size or smaller than the resolution errors. 
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However at 1/2° resolution, linear interpolation errors are noticeable, while cubic interpola- 
tion introduces errors which are small compared to the effect of resolution. 

4.b.5 Solutions with two observations 

A second “obs” at 45N, 50W, of 27 ms -1 from 350° is added to our nominal experiment. 
The two observations are similar to the two critical observations close to the center of the 
QE II storm studied by Hoffman (1982 [22]). The response of 2d-VAR to “a second obser- 
vation” (panel in Fig. 10) by itself is analogous to the results of the nominal case. That is, 
the response is aligned with the observed wind direction, and is proportional to the observed 
wind speed. Adding a second observation results in a solution due to “two observations” 
which is qualitatively in agreement with, but weaker than, the “sum of solutions” for each 
observation separately. However, the sum of the individual solutions is considerably weaker 
than the solution for the two observations. 

4.c Discussion 

Clearly there is room for improvements in the minimization procedures used. In other 
problems of this sort good preconditioning is vital for fast convergence. This may be helpful, 
but all variables in the present problem are wind components and have uniform scaling. 
Furthermore, the discussion regarding Fig. 5 shows that the number of iterations required 
must grow at least linearly with the number of grid intervals required to traverse the domain. 
Therefore the multigrid strategy of solving the same problem with increasing grid resolution 
initialized from the previous solution should be generally useful. A spectral representation 
might also be helpful in this regard. Alternatively, but for conventional observations only, a 
good initial estimate at any resolution might be obtained by a Green’s function approach. 
In this approach we would add approximations to the solution for each single “obs” to 
the background. For each “obs” a standard solution would be rotated to alignment with 
the “obs” and scaled by the difference between the “obs” and background to provide the 
approximate single “obs” solution. 
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The results described here have some bearing on a high resolution 2d-VAR. With low 
weight given to the background and high weight given to the filtering and dynamic constraints 
the response of 2d-VAR is spread out. Any resolution with a grid spacing small enough to 
resolve the response will give the same result when fully converged. A high resolution analysis 
requires greater data density and less weight given to the filtering and dynamic constraints. 
With this adjustment to the A weights, the response of 2d-VAR is spatially narrower. For 
this reason and to resolve the details potentially present in a high density data set, smaller 
grid spacing is required. 

5 Comparison to OI 

As mentioned earlier, 2d-VAR can be considered a form of OI in which the background 
error correlations are parameterized with a small number of degrees of freedom. Here we 
demonstrate that with properly chosen A weights, 2d-VAR produces analyses similar to those 
produced by OI. It should be noted that OI uses geostrophy as a strong constraint, while 
2d-VAR uses the time rate of change of vorticity computed from the full primitive equations 
(as described in § C) as a weak constraint. Thus differences might be expected in the tropics. 
The example shown here is therefore in the extratropics. 

Fig. 11 shows the NSCAT observations in the region used for the experiments described 
in this section. For the purposes of comparing assimilation techniques, the unique 25 km 
NSCAT winds are thinned to ~ 75 km resolution by retaining only every third WVC across 
track. Fig. 12 shows the background wind field used. NSCAT observation locations are 
marked by small black dots and the background wind speed is contoured every 2 ms F 

The implementation of OI used, the Theater Analysis Procedures (TAP) (Nehrkorn 2000, 
Nehrkorn and Hoffman 1996 [41, 42]), is similar to the Navy FNMOC OI of Goerss and 
Phoebus (1992 [16]). Standard values of the statistics were used here, treating the scat- 
terometer winds as ship reports with a default wind component error standard deviation 
(s u ) of 2.2 ms -1 . 
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Fig. 13 shows an 01 analysis for a region of rev 2607 at 2000 UTC 15 February 1997. For 
this analysis, s u has been increased by a factor of two to avoid overfitting the dense satellite 
observations. Fig. 14 shows the analysis increments (analysis-background) for the analysis 
in Fig. 13. The increments are significant and broaden the circulation around the cyclone as 
seen in the observations. Statistics for the background and analyzed wind field are found in 
Table 3. The rms vector analysis increment is 2.716 ms -1 . 

We then tried to match a 2d-VAR analysis to the OI solution for the same swath of data. 
In our search, we fixed the 2d-VAR A weights for background and smoothness constraints 
{{XvwmAlapAdivAvorAdyn) = (1,2, 8, 2, 8)), and varied the weight given to the data 
(Accwv = 5, 10, 20). 

Table 3 also shows the differences between the 2d-VAR and OI analysis presented above. 
A particularly good match, with rms vector wind difference of only 0.584 ms 1 was obtained 
using Accwv=5. 

6 Ambiguity removal approach 

6. a Definition of scatterometer “obs” functions 

For conventional ship or buoy observations J CONV is simply the sum of squared vector differ- 
ences between the “obs” and the analysis. For scatterometer data, we use two formulations. 
The first formulation, the method of H84, uses the retrieved scatterometer wind ambiguities 
in two “obs” functions — J AMB and J SPD . For a single location, when the analysis is close 
to one of the ambiguities, the misfit measured by J AMB approximates the squared vector 
differences between the ambiguity closest to the analysis and the analysis. Away from the 
ambiguities J AMB approaches a constant. The variability in wind speed among ambiguities 
at a single location is relatively small (H84). Therefore the observed mean or rms wind 
speed averaged over the ambiguities at each location is representative of the scatterometer 
wind speed. This value is used in the definition of J SPD . The second formulation uses the 
observed backscatter values in J NRC s- For a single location, the misfit measured by -Ames is 
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the likelihood function used to retrieve ambiguities. The likelihood function is equal to the 
squared normalized differences between the simulated and observed scatterometer a values. 
An optional log term may be included when the a 0 variance is parameterized in terms of the 
analyzed a 0 as explained below. Thepaut et al (1993 [56]) used a slightly simplified form of 

J NROS • 

Alternative formulations have been suggested for the loss function for scatterometer data. 
For example, operational use of ERS-1 and ERS-2 data at ECMWF uses the formulation of 
Stoffelen and Anderson (1997 [54]). Also, the variational QC of Andersson and Jarvinen(1999 
[1]) might be extended to scatterometer data, by treating each ambiguity as a ship report. 
These reports would be inconsistent at each location and the QC algorithm would essentially 
select one ambiguity at each location. 

Example plots of J NRCS , Jamb and the ECMWF formulation are given in Fig. 15 for 
a single WVC. The precise definitions of J N rcs, Ja mb an d Jspd are given in the following 
subsections. If a 0 values are used directly in 2d-VAR, the loss function due to a single WVC 
is the MLE (Fig. 15a). This surface shows the complexity of the highly nonlinear wind 
retrieval problem. Note that the center of the figure is not shaded because the values there 
are very large, O(10 5 ), compared to those plotted in the vicinity of the minima, 0(10°- 10 2 ). 
The MLE has multiple minima, four in this case, which correspond to the ambiguous wind 
solutions. The minima are distributed in an annulus of near-minimum values at roughly the 
same wind speed. The four minima are not equally deep, which implies that each solution 
has a slightly different likelihood of being the true wind. In fact, if there were no “obs” or 
model function errors, only one solution would be found the true wind. 

If ambiguous winds are used in 2d-VAR, then the loss function of a single WVC is Jamb 
(F ig. 15b). This approximation of J NRCS loses some of the structure seen in Fig. 15a, but 
still reflects the four wind solutions as minima and is less severely nonlinear, especially near 
the origin. Note that if the background wind field is close to one of the four wind solutions, 
the influence of other observations or meteorological constraints will be necessary to shift 
the analysis from near one minimum to another during the 2d-VAR minimization. 
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For comparison, we also show the two- wind loss fuction of StofFelen and Anderson used 
at ECMWF (Fig. 15c). This is the best behaved of the three functions, but is also the 
crudest approximation to the MLE. Only the two most likely winds are used since one of 
these solutions is very close to the true wind more than 90% of the time. The smooth nature 
of this function allows the analysis to be moved from one minima to another during the 
assimilation more easily than either the MLE or H84 formulation. However, this formulation 
may also permit an analysis solution which is not very close to either minima because of its 
smoothness. 

6.a.l Loss function for scatterometer backscatter values 

For each backscatter “obs”, the normalized departure is 

rr° — rr° 

(13) 

sa 

Here 

• <t° 0 is the backscatter observation. 

• a° a is the simulated backscatter observation. 

• S(j is the standard deviation of the a 0 “obs'’ . 

Theoretical estimates of the standard deviation of scatterometer measurements show that 

^ = K p a° (14) 

(Fischer 1972 [14]). Here K p is proportionality constant. The subscript p is used because K p 
is originally defined in terms of the radar power. The JPL model for the standard deviation 
of the NSC AT a 0 “obs” extends this model to 

si = K pA (a 0 ) 2 + K pB g° + K pC (15) 

(Long et al. 1988 [36]). (For ERS-1,2 (14) is used.) When evaluating s a , o° may be the 
observed value or the simulated value. If the simulated value is used, a logarithm term 
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should be added to the loss function to be consistent with maximum likelihood estimation 
theory. That is 

J NRCS T XX + lnM]. ( 16 ) 

The term in brackets comes about from taking the logarithm of minus the normal probability 
function of cj 0 a given < 7 ° 0 . In practice the logarithmic term is not normally used. 

A model function or forward model is used to simulate cr°. As described in §2, the 
model function calculates cr° as a function of wind speed and direction and the geometry 
and polarization of the observation. The geometry is specified in terms of the pointing 
direction of the antenna and the incidence angle. For NSC AT the NSCAT-2 model function 
is preferred. 

6. a. 2 Loss function for ambiguous scatterometer winds 

For ambiguous wind vectors the formulation of H84 is used. At each WVC, 

Jamb = £ SnSCAT II d l\ l ~ eX PM*/Jo)]> ( 17 ) 

k—\ 

where the sum is over all WVCs. and 

d 2 k = (u a - u k ) 2 + (fa — ffc) 2 - (16) 


Here 

• K is the number of ambiguities. 

• k is the index of the ambiguities. 

• (u k ,v k ) is the kth wind ambiguity. 

• d k is the magnitude of the vector difference between analysis and ambiguity k. 

• d 0 is the wind speed scale for the “obs” . 

do = Voh- ( 19 ) 
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Vo is the mean or rms wind speed for the “obs”. The mean or rms is over the K 
ambiguities, 


V„ = K-' ■£ <fM 2 + M 2 , 01 V 2 = K-' X>*) 2 + w 2 - (20) 

k = 1 fc=l 

• 7 is the scaling factor for the wind speed. 

• s N scat is the wind component standard deviation for the WVC. 

Nominal values for 7 and Snscat are 2 and 1 ms -1 respectively. 


6. a. 3 Loss function for scatterometer wind speed 

The wind speed operator is designed to take advantage of the relatively unambiguous wind 
speed information in scatterometer data. In analogy to (6) we define J SPD as, 

(V a ~ V 0 f 


dsPD — 


( 21 ) 




Here 


• 14 is the analyzed wind speed. 

• V 0 is the observed wind speed. 

• sk is the wind speed error standard deviation. 

This formulation has been used with SSM/I observations (Atlas et al. 1991, 1996 [2, 3]). 
The same loss function may be used for line-of-sight (LOS) wind speed in place of wind 
speed as long as the analyzed LOS wind speed can be calculated. This requires knowledge 
of the antenna pointing direction. In H84 a different wind speed functional was introduced. 
It is equivalent to (21) if sy is multiplied by the factor ^1 + V 0 f\ 0, where V 0 is a constant 
2 ms -1 . Note that this has the effect of giving more weight to low wind speed reports. 
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6.b Scatterometer quality control 
6.b.l Gross checks 

Gross checking is performed when the data is read from source files and prevents nonsensical 
or nonphysical values from being stored in 2d-VAR data structures. For example, NSCAT 
values of < 7 ° > 30 dB are eliminated. Small negative values of a 0 are nonphysical but are 
indicative of very light winds (Pierson 1989 [46]), and in 2d-VAR are replaced with a small 
positive value equivalent to —60 dB. 

6.b.2 Background and analysis checks 

Background and analysis checks are done by comparing observations with the current analy- 
sis. At the start of the process, when the current analysis is the background, this procedure 
is a background check. Later it is a preliminary analysis check. Data failing the background 
check may be considered by a later (and stricter) analysis check. This allows some of the 
data eliminated by earlier QC back into the analysis. 

As in H84, if the magnitude of the difference between wind vectors (or wind speeds) is 
greater than the average of the observed and analyzed wind speed, the observation fails QC. 
If the difference between the wind directions is > 60°, the observation fails QC. However, 
very small or zero winds are handled as a special case, so that when both observation and 
analysis indicate calm conditions, the previous checks are bypassed. The background and 
analysis checks apply to conventional data as well as to scatterometer winds. For ambiguities, 
only the ambiguity closest (i.e. , with the smallest vector difference) to the current analysis 
is checked. 

For NSCAT ct° data, we use a simple departure test for each cr° measurement in dB space. 

<J°o - O°o. > 7 NSCAT ( 22 ) 

We have used Jnscat — 9 dB for both background and analysis checks, which typically 
eliminates ~ 5% of the data. 
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7 Use of ambiguous winds versus a 0 observations 

When using a 0 observations, 2d-VAR essentially carries out a simultaneous retrieval of 
NSCAT winds for all data presented to the analysis. If NS CAT ambiguous winds and 
< 7 ° data represent the same geophysical information, namely ocean surface wind stress, w'e 
expect 2d-VAR to produce similar wind analyses using either data source. If this is correct, 
A for ambiguous wind and < 7 ° data must exist which produce similar analyses. To determine 
such equitable A for ambiguous winds and a 0 data, we can examine any measure of the 
misfit between data and analysis versus a range of A. We say the A are equitable when the 
misfit of the analysis of the ambiguous winds is the same as the misfit of the analysis of the 
( 7 ° data. This provides an example of tuning the A weights. Note that differences between 
winds-based and cr°-based analyses will also arise for other reasons: 

1. Quality control decisions made for cr° and wind data may result in different data being 
used, 

2. WVCs with no ambiguities may contain valid cr° data, and 

3. The ambiguous wind observation loss function, J A MB , is an approximation to the cr° 
data loss function (J N rcs) (compare Figs. 15a and 15b). 

The first two reasons will only effect a small fraction of the data presented to 2d-VAR while 
the third reason may effect all data to a small degree. With respect to the second reason, 
note that a single a 0 value is not sufficient to infer a set of ambiguities but is sufficient to 
refine the background wind in 2d-VAR. With respect to the third reason, the a 0 loss function 
(J NRCS ) will permit a wider variety of wind directions in the minimizing analysis compared 
to the ambiguous winds loss function ( J AM b)> because of the shapes of the functions surfaces. 
That is, J NRCS has an annulus of minimum values (all wind directions) with secondary local 
minima embedded while the minima in J AMB are more clearly separated. 

To evaluate the analysis fit to the data, we use the rms wind speed difference between 
the 2d-VAR analysis interpolated to NSCAT WVCs and the first ranked (i.e., highest MLE) 
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JPL ambiguity. For two North Atlantic cases, one on 19 September 1996 of Hurricane Lili in 
the western Atlantic and one on 27 October 1996 of a cyclone west of Ireland, analyses using 
ambiguous winds and a 0 data were created using a range of A. In these experiments the 
nominal setup includes a 1 x 1° grid over an area 20° of latitude by 20° longitude centered 
over an NSC AT data swath. The analyzed field is interpolated to observation locations by 
linear interpolation. The A for background constraints are held fixed for all analyses at 


(AvWM) A lap 7 A D iv, AvOR? A DY n) 


(1,1,4,1,16)^. 


The minimizations were allowed to run until the gradient test was satisfied (e — 10 3 ). 1 he 
A weights used for the winds and a 0 analyses are 


A A mb,nrcs - 1/1024, 1/256, 1/64, 1/16, 1/4, 1, 4, 16, 64, 256, 1024. 

Smaller A weights for scattterometer data required fewer iterations to reach convergence 
(typically 10-30), and larger A weights required more (200-600). The stronger constraint 
imposed by larger A weights (> 4) often caused the minimization to fail, since a new search 
direction could not be found to reduce the cost function any further. 

Fig. 16 shows the rms analysis-NSCAT observations difference as a function of A weight 
for both ambiguous wind (w) and <r° (s) analyses for the October 27 case. For very small 
A, the analyses of wind and < 7 ° data are scarcely changed and fit the wind observations 
equally poorly (~ 4 ms~ l ). For large A (> 4), the analyses of wind and a 0 data fit the wind 
observations much more closely, and again, equally (~ 1.5 ms* 1 ). The lower limit of the 
rms fit is governed by our choice of A 2 for the background constraints. If the background 
constraints were weakened (i.e., smaller A VW m> A lap , A div , A VO r)> the fit of the analysis to 
wind observations would improve. In the limit of no background constraints, the analysis 
would fit the observations exactly. But without any constraints on smoothness or divergence 
and vorticity, the solution may be unrealistic. The lower limit of rms fit is, in this case, 
also dependent on the robustness of the minimizer. The minimizations in analyses with 
A nrcs> 1/4 and A amb> 256 fail to converge sufficiently because a direction of descent 
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cannot be found to within machine precision. This explains why the fit to a 0 does not 
improve for A nrcs > 1* 

A background field, temporally interpolated to 2230 UTC October 27, 1996 from 6- 
hourly NCEP 1000 hPa wind analyses, was used for the tuning experiments presented next 
(see Fig. 17). (The other tuning case - Hurricane Lili, October 17, 1996 - produced similar 
results and will not be shown.) To illustrate that ambiguous winds and a 0 data represent 
the same information, Fig. 18 shows analyses and increments for \amb = 4 and X^ RCS =lf4, 
The upper panels show 2d-VAR analyses and the lower panels show analysis increments 
times ten. The analyses are very similar in pattern and magnitude. While the analysis 
domain rms vector magnitude is 3.144 ms~ l for the winds-based analysis increments and 
2.198 ms~ l for the (j°-based analysis increments, it is 1.856 ms~ l for the difference of the 
two analyses. Differences between the analyses are subtle. The center of the low is displaced 
very slightly to the east in the ambiguous winds analysis relative to the background and 
the a 0 analysis. Also, in the region south of Ireland, winds from the southwest in the 
background are changed to more southerly winds in the ambiguous winds analysis (easily 
seen in the increments, lower panels). This change produces a circulation pattern around 
the cyclone which is slightly less circular than in the a° analysis. This is very likely due the 
the difference in the observation loss functions. J A MB draws the solution toward one of the 
wind ambiguities, generally not permitting a solution to fall between minima. In contrast, 
T N rcs tends to allow any solution within the annulus of minimum values, with only a minor 
preference for the retrieved solutions (i.e., for the local minima within the annulus). Apart 
from these minor differences, it is clear that wind and a 0 data, while handled very differently 
by 2d-VAR, produce nearly the same result given an appropriate choice of A amb and A nrcs- 

8 Dual ambiguity processing 

While more than two winds are often retrieved from a group of colocated a 0 measurements, 
one of the two ambiguities with the highest MLE values is very likely (>90%) to be closest 


29 



to the true wind (JPL 1997 [32]). This implies that <10% of the third and fourth ranked 
ambiguities are closest to the true wind. This information can be used to QC the data given 
to 2d-VAR, and improve ambiguity removal. 

We observed that using all available ambiguities in 2d-VAR, or equivalently using the 
< 7 ° values, leads to patches of poor ambiguity selection. This is especially true where there 
are four ambiguities and the direction of the the third or fourth ambiguity is close to the 
direction of the background wind field. Using “dual ambiguities” (i.e., the two ambiguities 
with the highest MLE values) narrows the choice to one of the two most likely possibilities. 
This allows for larger (> 45°) wind direction corrections to the background field. Using dual 
ambiguities only, however, naturally limits ambiguity selection to the first or second ranked 
ambiguity since the analysis will generally be drawn toward one or the other. 

To allow for selection from all ambiguities, we have developed a two stage analysis pro- 
cedure which blends the use of dual and all ambiguities. In the first stage, we use dual 
ambiguities for 50 iterations of the minimizer to draw the analysis toward one of the two 
most likely observed winds. In the second stage, we include all ambiguities and continue to 
minimize until the convergence criteria are met. An alternative is to switch to a 0 values after 
the first stage. The second stage allows for readjustment of the analysis in the few locations 
where the third or fourth ranked ambiguity is closest to the true wind. 

During the first stage, we only use dual ambiguities which pass dual ambiguity quality 
control or “dual QC”. Dual QC is based on a conceptual model of dual ambiguities. Nom- 
inally, we expect the wind directions of dual ambiguities to be opposed by ~180.° If not, 
the retrieved winds are suspect and are not included in the first (dual ambiguities) stage. 
Specifically, retrieved winds fail dual QC and are not used in the first minimization if the di- 
rections of the two most likely ambiguities are within 37 t/ 4 of each other. Dual QC typically 
eliminates ~20% of WVC’s from the first minimization. 
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8. a Dual QC ambiguity removal examples 

One measure of the performance of ambiguity removal with 2d-VAR is a direct comparison 
with another successful ambiguity removal scheme. Here we compare our results with the 
NWP-initialized or “nudged” median filter used by JPL for ambiguity removal, and with 
ambiguities chosen to be closest to the NWP analysis used as the background. 2d-VAR 
winds and median filter winds agree in a large proportion of the WVCs, typically 95%. 
Detailed statistical comparisons for the whole NSCAT mission are given by Henderson et al. 
(2001 [ 21 ]). 

In cases of disagreement both sets of winds are still horizontally coherent. In our synoptic 
evaluation, we find 2d-VAR winds tend to be more meteorologically reasonable. Recall that 
both methods are nonlinear and have multiple solutions. Because 2d-VAR requires dynamic 
consistency and the median filter does not, solutions of the median filter are not necessarily 
solutions of 2d-VAR. For example, two blocks of winds with opposite directions can satisfy 
the median filter, but not 2d-VAR. 

As a test of our dual-QC procedure we initially processed 234 orbits of NSCAT data based 
on the NSCAT-1 model function during the period 15 October-5 November 1996. JPL and 
2d-VAR ambiguity selection differences are predominately at lower wind speeds (Table 4). 
To further assess the performance of 2d-VAR, we use GOES-8 imagery to locate mesoscale 
meteorological features in the North Atlantic, for the subset of 29 orbits passing over the 
North Atlantic when the sun elevation angle was sufficient to make use of the GOES visible 
channel (Grassotti et al. 1999 [20]). We examined 11 North Atlantic scenes where JPL and 
2d-VAR ambiguity selection differed. From the 11 scenes, 29 subregions were examined in 
detail. For brevity, we present one of the subregions below and then summarize the results 
of all 29. 

8.a.l Case of 26 October 1996 

Fig. 19 shows GOES-8 IR brightness temperature over the North Atlantic at 1315 UTC 
26 October 1996. NSCAT data locations are marked by the fine black dots, and WVC’s 
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where JPL and 2d-VAR-selected ambiguities differ are marked by heavy blue dots. The 
scatterometer crossed 30N at 1314 UTC. Differences between JPL and 2d-VAR-selected am- 
biguities occur in small patches near (35W, 29N), (37W, 41N), (30\V, 45N), and (32W, 57N) 
in Fig. 19. In region A, satellite imagery shows that the scatterometer passed over a mature 
cyclone at 36W, 41N. 

Fig. 20 shows ECMWF (thin red barbs) and 2d-VAR (thick black barbs) analyses in 
region A over an image of GOES-8 IR brightness temperature. NSCAT WVC locations are 
marked by small black circles, and points where JPL and 2d-VAR ambiguity selection differ 
are marked by larger blue circles with white centers. It is interesting to note several effects 
of the NSCAT data on 2d-VAR surface wind analysis. The circulation center in 2d-VAR 
analysis is moved north and east of the center in the ECMWF analysis. The placement in 
the 2d-VAR analysis is in better agreement with the satellite image. Also, the wind speeds 
around the cyclone are higher (roughly doubled) in the 2d-VAR analysis than in the ECMWF 
analysis. 

All JPL ambiguities in subregion A1 of Fig. 20 are plotted in Fig. 21. The ambiguous 
winds are consistent with a cyclone centered at 36W, 41N. Fig. 22 shows JPL and 2d-VAR 
selected ambiguities for subregion Al. A physically unlikely east- west wind shift exists in 
the JPL ambiguites in the northwest quadrant of the storm (37W, 41. 5N). The satellite 
imagery is more consistent with the circular wind flow around the cyclone center evident in 
Fig. 22b, which shows 2d-VAR selected ambiguities. Note that 2d-VAR positions the center 
of circulation > 50 km south of the wind shift in the JPL selected ambiguities. 

8.a.2 Summary of All Cases with Colocated GOES Imagery 

The subregion presented in §8.a.l is one of 29 cases where 2d-VAR and JPL ambiguity 
selection differ that were examined for the three week period. These cases include cloudy 
and clear scenes, tropics and mid-latitudes, and coastal and open water regions. The re- 
sults from each subregion are categorized by the differences between the JPL and 2d-VAR 
selected ambiguities. Using GOES-8 imagery and general knowledge of the local wind field 
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from ECMWF analyses, 2d-VAR selected ambiguities are determined to be either 1) clearly 
improved over JPL selected ambiguities, 2) different from JPL but equally plausible, or 3) 
clearly worse than JPL selected ambiguities. For reference, 2d-VAR selected ambiguities 
were determined to be a clear improvement over JPL selected ambiguities for the subregion 
presented in §8.a.l. For the 29 subregions we examined, 

• 10 subregions are clearly improved over JPL selected ambiguities 

• 16 subregions are different from JPL but equally plausible 

• 3 subregions are clearly worse than JPL selected ambiguities 

Cases classified as “equally plausible” often occur when the evidence from the imagery is 
not decisive. For example, sometimes a front is aligned with a cloud feature in both 2d-\AR 
and JPL selections but displaced one or two WVCs in the perpendicular direction. In other 
“equally plausible” cases all ambiguities are inconsistent with the imagery. 


9 Concluding remarks 

2d-VAR provides an alternative method of ambiguity removal for NSCAT data. In this paper 
we describe both 2d-VAR and its application to NSCAT data. Experiments with a single 
ship observation illustrate the numerical behavior and sensitivity of 2d-VAR. Note that single 
observation experiments with 2d-VAR when normalized to have a unit response at the “obs” 
location, are the background error correlations with the single “obs”. (This is a general 
result: Thepaut et al. (1993 [56]) essentially showed this for a more complicated situation.) 
In experiments based on no tuning other than the subjective tuning described by H84, the 
response of 2d-VAR to a single ship observation is very similar to the response of 01 analyses 
making use of models of forecast error statistics used operationally by FNMOC (Goerss and 
Phoebus 1992 [16]). Therefore the subjectively tuned 2d-VAR approximates 01. Idealized 
01 experiments with a swath of unique NSCAT winds demonstrate the close relationship 
between 2d-VAR and 01. 2d-VAR analyses can be used to select ambiguities by choosing 
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the ambiguous wind closest to the analyzed wind field at each WVC location. Comparisons 
between the use of ambiguities and backscatter in 2d-VAR show the near equivalence of 
these two data for this purpose and also demonstrate tuning the A weights. Whether using 
ambiguities or backscatter, the objective function is highly nonlinear — multiple solutions 
are possible. The analysis tends to settle on the ambiguity closest to the background field. 
When the background quality is not sufficient, the fraction of third and fourth ambiguities 
chosen is larger than expected: NSCAT and buoy colocations show that more than 90% of 
selected ambiguities should be the first or second ambiguity. This leads to the concept of 
“dual-ambiguities” quality control (dual-QC). 

Ambiguity removal using 2d-VAR show's promise and may have advantages over median 
filter techniques. While JPL and 2d-VAR ambiguity selection differ for only ~5% of WVC’s, 
the differences tend to occur in patches. Patches of poorly selected winds can be more harmful 
to analysis and forecasting systems than scattered errors if not properly quality controlled. In 
~30% of the cases examined, 2d-VAR selected ambiguities are a clear improvement over JPL 
selected ambiguities. In ~60% of these cases, there was either insufficient information to 
determine which choice was more likely or a reasonable solution w T as not present in the 
ambiguities. In the remaining ~ 10% of these cases, the JPL selected ambiguities are 
preferred. In a companion article (Henderson et al. 2001 [21]), w r e process the entire 9- 
month mission using the latest version of the data from JPL based on the NSCAT-2 model 
function. 

Further work could be done to investigate the use of the scatterometer backscatter mea- 
surements (ct°) directly for ambiguity removal. This may prove most useful in cases when 
a suitable solution does not seem to exist in the retrieved winds. Many of the principles 
learned here are applicable to NASA’s current Seawinds scatterometer on the QuikSCAT 
satellite. In addition there are other factors which we have studied in other contexts which 
impact ambiguity removal. Positional errors in the background field can have a substan- 
tial effect on ambiguity selection. Henderson et al. (2001 [21]) show the effect of using a 
very poor background. Smaller positional errors might be handled wdthin 2d-VAR using our 
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method for feature calibration and alignment (Hoffman et al. 1995, Hoffman and Grassotti 
1996, Grassotti et al. 1999 [25, 24, 19]). Rain effects seem to be important for QuikSCAT 
(Jones et al. 2000, Mears et al. 2000 [31, 37]). However for NSC AT, limited experiments 
with rain flagging show small impacts on 2d-VAR results. Grassotti et al. (1999 [20]) de- 
veloped a rain flag for NSCAT data based on geostationary visible-infrared imagery. That 
study was a demonstration restricted to three weeks of data in the fall of 1996 in the North 
Atlantic. For most cases studied by Grassotti et al., the effect of removing flagged WVC on 
ambiguity removal using 2d-VAR was small. Still, occasionally, large differences can result 
from eliminating just one WVC. 
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A Representation of the wind fields 

For the purpose of minimization, we consider J to be a function of X, which is known as 
the control variable. We must be able to map the control variables to a complete gridded 
wind field. That is 

y a = A4(X). (23) 

Here 

• V a is the gridded wind analysis. 

• A4 is the mapping operator from control variables to gridded winds. 
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X is the control variable vector. 


Nominally V a is the 10 m neutral stability wind field. All observed winds should be translated 
into a 10 m neutral stability wind. If only the height of the observation is known, then the 
neutral logarithmic wind profile may be used to adjust the wind observation to 10 m. The 
control variables might contain grid point wind components, spectral coefficients of anal} sis 
increments, or some other representation of the wind fields. In one sense X must completely 
describe the wind field, however any constant auxiliary information may be used by M. 

The grid is specified in terms of longitude and latitude denoted for convenience by the 
variables x and y, respectively. The usual setup is to evaluate all the functionals over a 
rectangular region in x and y. However the data acceptance window for evaluating the J 0 
terms and the integration domain for evaluating the Jb terms may be specified individu- 
ally. (In fact one could specify these domains individually for each component of these cost 
functions.) To begin, we define the integration domain: 

For both x and y there is a starting location, an increment S, and a number of grid 
boxes tl. These grid boxes are numbered from 1 to n and the corresponding grid points are 
numbered from 0 to n, 


Grid box 
Grid point 


1 


2 


n 


0 1 2 n - 1 n 


To evaluate some of the finite difference operators and to allow higher order horizontal 
interpolation, a boundary zone one grid box wide is added to the grid. The grid point 
indices therefore run from -1 to n + 1, and the boundary grid boxes are numbered 0 and 
n - f 1 . Therefore 


Grid box 
Grid point 


0 


n 


n + 1 


-1 


0 


1 


n — 1 


n 


n+ 1 


Grid points 0 and n at the integration boundary may be either active or passive, i.e. allowed 
to vary during the minimization or held fixed. If they are held fixed they are part of the 
boundary. In the experiments reported here they are active and part of the control vector. 
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The domain may be periodic in the x direction. In this case grid points i and i + 2i x/8 are 
the same. Our implementation does not allow pole points. 

With the grid specified, we now describe the operator M. The fixed points of the grid 
are held constant or reset to their fixed values. All other values are copied from locations 
in X to locations in V^. Therefore X may be taken to be the concatenation of portions of 
the gridded wind component arrays which have been reshaped into vectors. The relevant 
portions are free to vary and do not extend beyond one period in the x direction. If the grid 
is periodic some values may be copied twice. 


B Horizontal interpolation 


The grid point values of the (u, v ) wind components are interpolated bilinearly to the obs 
locations. (That is, u and v are separately interpolated first in longitude and then in lati- 
tude.) In some experiments, a bicubic interpolator is used. Since the interpolated value is 
a linear function of the data values at the grid points, we simplify the adjoint calculation 
by formulating the interpolation as a weighted sum of grid point values. The key point 
is that these weights are constants depending only on the “obs ’ location. Interpolating in 
one dimension and then in the other, shows the weight for a grid value is the product of 


the weights for interpolating in each single dimension. In a single dimension, the weights 
are most easily determined from the Lagrange multiplier representation of the colocating 
polynomial, 


/ 0 ) = 


n 


X - Xi 


it) x o - x i 


fv 


(24) 


where ft are the data at locations x y The sum and product in (24) are from 1 to N + 1, 
where N is the order of the polynomial. Cubic interpolation for example requires four data 
points. When x is equal to one of the x y the weights in (24) reduce to one for j and zero 
otherwise. The bicubic interpolation actually uses a 12 point mask, following Ritchie et al. 


(1995, Fig. 2 [48]). 
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C Calculation of the background constraints 

The background constraints are all formulated as 

Jm = 4~ f (Km(V a ) - n m {V b )fdA. (25) 

That is, each background constraint term, J m , listed in Table 1 may be written as the scaled 
integral of the squared difference of some operator applied to the analysis field and to the 
background field. 

Below, several of the operators 7Z m are put in the form, 

n = v-Q. ( 26 ) 

In numerical experiments, we found that this version of the Laplacian leads to an instability 
in the minimization. Instead the Laplacian is calculated at each grid point and then J VW m 
and J LAP are calculated in the same way, by averaging four adjacent grid point values to 
estimate the value at the center of the grid box. (See the discussion in § C.2.) 

The filtering constraints are based on minimizing the vorticity, divergence and Laplacian 
of the analysis increments. We call these constraints filtering because they tend to eliminate 
small spatial scales. The dynamic constraint is based on minimizing the analysis increments 
of the time rate of change of vorticity. An second similar dynamic constraint could easily be 
formulated in terms of divergence, but would require introducing surface pressure into the 
problem. The background constraints and scales are specified in terms of S m , 7Z m and Q m 
as follows: 


m 


dim 

Qm 

VWM 

t~ 2 l 4 

V 

— 

LAP 

T ~ 2 

(V 2 u, V 2 t») 

— 

DIV 

T-2 L 2 

<1 

V 

VOR 

T~2 L 2 

-V • k X V 

k x V 

DYN 

T~ 4 L 2 

dQ/dt 

X 

l 

+ 

T 
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Here 


• / is the Coriolis parameter. 

• F is the surface friction term in the momentum equation. 

The definitions of the Q m for J DIV and J VOR follow immediately from the definitions of 
the Tl m and the vector identities. To compute the dynamic constraint and to determine the 
<3 dyn , we start with the equations from Bourke et al. (1977 [6]) for the time rate of change 
of vorticity. We then eliminate horizontal (but not surface) friction and vertical advection 
terms. Horizontal friction is negligible at synoptic scales. Formally the vertical velocity in 
sigma coordinates is zero at the surface. 


C.l Finite differencing 


Simple finite 
geometry, 


difference forms are used. If the components of Q are (p, q), 

tan $ 

V • Q « p x + Q y — 


then in spherical 


(27) 


and 


V Q — (?i) Qy ) > 


(28) 


where for any p and q, 

_ P i±hj ~ Pi-hj n ■ ■ - ftd+i ~ (29) 

2a cos yjSx ’ 2a8y 

These finite difference forms are applied directly at a grid point to evaluate <. However, the 
71 are evaluated at the center of the grid box. Instead of applying (29) directly in (27) at 
the grid points and then averaging, we apply (27) at the center of the grid box. We evaluate 
% by first averaging in y and then differencing in x. The overbar indicates a value at the 
center of the gridbox. Similarly, we evaluate q y by first averaging in x and then differencing 
in y. For q, we again take the average of the corner grid point values. 
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C.2 Evaluation of integrals 

All integrals are put in the form of (25). For J v WM an d th ere are tw0 integrals of this 
form— one for each wind component. 

For all J m we discretize the integral as 

j m = < 3 °) 

where 

• A Rij is the estimate of R m {V a ) - TZ m {V b ) at the middle of grid box ij. 

• Aij is the area of grid box ij. 

As mentioned for J VWM and J LAP , the Rn are estimated as the average of the four corner 
grid point values. For the other constraints, the R m are estimated as described in §C.l. 

C.3 The gradient calculation 

Simple finite difference formula are used because J is ad hoc and the A are adjustable. Thus 
we may easily subsume the approximation errors introduced by the finite difference and 
integration formula into the assumption that minimizing J provides the optimal analysis. In 
other words, we identify the best analysis as the minimizer of the discrete formulation of J. 
These considerations suggest choosing simple averaging and difference forms. However, it is 
vital that once the finite difference version of J is chosen, the gradient of J be calculated as 
exactly and precisely as possible. 

The calculation of J b and its gradient may be made efficient, both in terms of calculations 
and memory, as follows. Since both the analysis and background fields are gridded, it is 
possible to evaluate all of the constraints by sweeping from south to north in latitude. In 
addition the total background function may be considered to be the weighted sum of many 
partial cost functions, one for each constraint and each grid box. Since the weights are 
fixed and known a priori, it is possible to accumulate the contributions of these partial 


40 


cost functions to the sensitivities of the total background function in the same sweep over 
latitude. 
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TABLE LEGENDS 


Table 1: 


Table 2: 


Table 3: 


Table 4: 


Summary of “obs” functions and background constraints used in the variational anal- 
ysis. Subscript a refers to analyzed values (e.g., V a is the analyzed wind vector), 
subscript b to background values, subscript o to observed values, and subscript k to 
the kth ambiguity retrieved from scatterometer data. Other subscripts are defined in 
the table, and x, and C are velocity potential, streamfunction and relative vortic- 
ity , respectively. For clarity scaling factors are not included in this table. Complete 
definitions and our implementation are detailed in the text. 


Summary of the A weights and results for the sensitivity experiments described in 
the text. In all experiments listed A VWM = 1 -In addition to the values of the other 
A weights the table lists the total number of function evaluations (N f ), the rms 
gradient of the objective function (V J(xl0 6 )), and the objective function itself (J). 
The analysis speed and direction (9) are given at the obs location. 

Statistics comparing surface wind analyses from optimal interpolation (01) and 2d- 
VAR with the background and the first analysis in the table. The 01 analysis uses 
twice the nominal wind component error standard deviation (s u ). In the remaining 
analyses, the 2d-VAR A weights are held fixed, excepting Xconv which is varied 
systematically. The statistics calculated are the mean and the rms of the wind speed 
for the held, the rms of the magnitude of the vector differences between the analyses 
and background (denoted increment^, and the rms of the magnitude of the vector 
differences between the analyses and 01 analysis (denoted difference^. 

Differences between JPL and 2d-VAR ambiguity selection as a function of wind speed 
for 15 October-5 November 1996. The column labeled N gives the number of com- 
parisons in each wind speed bin. 
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FIGURE LEGENDS 


Fig. 1: Analyzed wind field for 1 /2° resolution, fully converged for the nominal case described 

in the text. Here, and in the figures which follow, the wind barbs are in knots and all 
results are presented on a 1° grid. Also the color scale shown here fin ms ) is used 
in the figures which follow. 

Fig. 2: Statistical interpolation results using NOGAPS statistics. 

Fig. 3: The effect of resolution. Results for 2, 1, 1/2 and 1/4° resolution. The panel labeled 

“1/2 degree resolution” repeats the results shown in Fig. 1. 

Fig. 4: Differences, scaled by 10, from the result for 1/8° resolution, for the results for 2, 1, 

1/2 and 1/4° resolution (shown in Fig. 3). (In the first panel, and in future figures, 
results greater than 20 ms 1 are colored (yellow) as if they were 20 ms .) 

Fig. 5: The sequence of results for increasing number of iterations for 1/2° resolution, starting 

from a zero initial estimate. 

Fig. 6: Differences (times 10) between the result using only 25 iterations at each resolution 

and the “fully” converged solutions ( shown in Fig. 3) for 2, 1, 1/2 and 1/4 resolution. 

Fig. 7: Solutions for different A weights. The nominal case is shown in the lower left panel. 

The X LAP , X DIV , A V0R , and X DYN increase by a factor of 4 in the right column, while 
A CCW v increases by a factor of 4 in the upper row. 

Fig. 8: Solutions for setting some of the X weights to zero. The nominal case is shown in the 

lower left panel. The X DYN is set to zero in the right column, while X Drv and Xvor are 
set to zero in the upper row. 
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Fig. 9: Differences ( times 10) from the solution for 1/4° resolution for runs at 1° and 1/2 

using bilinear and bicubic interpolation procedures. 

Fig. 10: Effect of adding a second surface wind observation. The difference ( times 10) of the 

solution for two observations minus the sum of the individual solutions is shown in 
the lower right panel. 

Fig. 11: JPL-selected ambiguities from rev 2607. 

Fig. 12: Background wind field, time interpolated to 1930 ETC February 15, 1997 from six- 

hourly ECMWF 10 m wind fields. Wind speed is contoured every 2 ms -1 and NSC AT 
“obs” locations are plotted as small black dots. 

Fig. 13: Wind field analysis from TAP with the standard deviation of “obs” error inflated by 

a factor of two. Wind speed is contoured every 2 ms -1 and NSCAT “obs” locations 
are plotted as small black dots. 

Fig. 14: Wind field analysis increments from TAP corresponding to Fig. 13. Wind speed is 

contoured every 2 ms~ l and NSCAT “obs” locations are plotted as small black dots. 

Fig. 15: 2d-VAR observation operators evaluated for one WVC using o° values directly (upper 

panel), and ambiguous winds (middle panel). For reference, the corresponding two- 
wind observation operator in use at ECMWF is shown ( lower panel ) (Stoffelen and 
Anderson (1997 [54])). 

Fig. 16: The rms fit of analyses to NSCAT wind observations as a function oflog 10 (X) weight 

for the October 27, 1996 case (rev 1025). Curves for analyses using ambiguous wind 
(w) and o° (s) observations are shown. 
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Fig. 17: NCEP background for winds and o° analyses, valid at 2230 UTC October 27, 1996. 

Fig. 18: Comparison of VAM analyses from wind and o° observations for A amb—4 and 

A nrcs=1/4, valid at 2230 UTC October 27, 1996. The upper panels are analyses 
and lower panels are analysis increments ( analysis-background times ten). 

Fig. 19: North Atlantic overview at 1315 UTC 26 Oct 1996 including GOES-8 IR brightness 

temperature (K). NSCAT WVC locations are marked by fine black dots. Heavy blue 
dots mark WVC’s where JPL and 2d-VAR ambiguity selection differ. 

Fig. 20: ECMWF ( thin red barbs) and 2d-VAR (thick black barbs) analyses in region A with 

GOES-8 IR brightness temperature. See Fig. 19 for the location of region A. NSCAT 
WVC locations are marked by fine black dots. Blue circles with white centers mark 
WVCs where JPL and 2d-VAR ambiguity selection differ. 

Fig. 21: All JPL ambiguities and GOES-8 IR brightness temperature in subregion Al. See 

Fig. 20 for location of subregion Al. Blue circles with white centers mark WVCs 
where JPL and 2d-VAR ambiguity selection differ. 

Fig. 22: JPL (a) and 2d-VAR (b) selected ambiguities, and GOES-8 IR brightness temperature 

in subregion Al. See Fig. 20 for location of subregion Al. Blue circles with white 
centers mark WVCs where JPL and 2d-VAR ambiguity selection differ. 
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TABLES 


Table 1: 


Table 2: 



Term 

Expression 

Description 


“Obs” 

functions 

CONV 

£(V a - V 0 f 

for the conventional data 

J AMB 

E rifc{l ~ exp[— CE, — Vfc) 2 ]} 

for the ambiguous winds 

JsPD 

£[|V„| - rms(|V 0 |)] 2 

for the scatterometer wind speed 

^NRCS 

£(A - A) 2 

for the backscatter 

Background constraints 

J VWM 

3 

i 

to 

on the vector wind magnitude 

*^LAP 

/[V 2 K-u 6 )] 2 + /[V 2 (u q -u 6 )] 2 

on the Laplacian of the wind components 

J DIV 

/[V 2 (Xa - Xb)] 2 

on the divergence 

JvOK 

/[V 2 ^-^)] 2 

on the vorticity 

J DYN 

KdCa/dt - dCb/dtY 

on the vorticity tendency 


Run 

^CONV 

^LAP 

^DIV 

^VOR 

^DYN 

\V\ 

e 

N f 

YJ 

J 

Nominal 

80 

4 

16 

4 

64 

18.65 

210 

994 

127 

12937 

“Obs” + 

320 

4 

16 

4 

64 

25.73 

210 

1060 

1822 

18715 

Constraints + 

80 

16 

64 

16 

256 

11.78 

209 

1060 

378 

21398 

Background — 

320 

16 

64 

16 

256 

20.83 

210 

1060 

39460 

41044 

w/o DIV, VOR 

80 

4 

0 

0 

64 

21.91 

210 

1060 

2188 

9278 

w/o DYN 

80 

4 

16 

4 

0 

21.64 

210 

498 

119 

10036 

LAP only 

80 

4 

0 

0 

0 

25.17 

210 

513 

107 

5793 
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Table 3: 


Name 


mean 

rms increment difference 

Background 


8.788 

9.287 

— 

— 

01 analysis, s u x 2 


9.508 

10.187 

2.716 

— 

2d-VAR analysis, A conv 

=5 

9.510 

10.183 

2.726 

0.584 

2d-VAR analysis, A conv 

=10 

9.601 

10.277 

2.942 

0.739 

2d-VAR analysis, A conv 

=20 

9.675 

10.351 

3.109 

0.935 


Table 4: 


Wind Speed (ms *) 

N 

% Different 

0-2 

260521 

18.379 

2-4 

836017 

11.169 

4-16 

5867770 

3.807 

>16 

142815 

1.438 

All 

7107123 

5.141 
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Fig. 22: 
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